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Abstract 

In the present work, we address a class of Cahn-Hilliard equations characterized by a singular 
diffusion term. The problem is a simpified version with constant mobility of the Cahn-Hilliard- 
de Gennes model of phase separation in binary, incompressible, isothermal mixtures of polymer 
molecules. It is proved that for any final time T, the problem admits a unique energy type weak 
, solution, defined over (0, T). For any r > such solution is classical in the sense of belonging to 

a suitable Holder class over (r, T), and enjoys the property of being separated from the singular 
^*) ' values corresponding to pure phases. 
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1 Introduction 

This paper is devoted to the mathematical analysis of the following class of parabolic systems: 

u t -Aw = 0, (1.1) 
to = -a(u)Au - ^-t^|Vu| 2 + /(it) -\u + eu t , (1.2) 

on (0, T) X O, fi being a bounded smooth subset of M 3 and T > an assigned final time. The system 
is coupled with the initial and boundary conditions 

tf| t= o = uo, in 0, (1.3) 
d n u = d n w = 0, on dtt, for t G (0,T) (1.4) 

and represents a variant of the Cahn-Hilliard model for phase separation in binary materials. Here, 
A and e are nonnegative parameters, where the case e > accounts for a viscosity effect that may 
appear the frame of Cahn-Hilliard models (see, e.g., [12]). Moreover, the function / stands for the 
derivative of the so-called logarithmic potential 

F(r) = (l-r)log(l-r) + (l + r)log(l + r), re [-1,1], (1.5) 
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namely, 



/(r) = log(l + r)-log(l-r)=logi±I, re (-1,1). (1.6) 



Finally, we assume that the function a has the form 

<r) = = f'(r). (1.7) 

1 — r z 

The initial-boundary value problem given by (|l.ip - (|1.3p will be noted as Problem (P) in the sequel. 

The above problem is a simplified version (assuming constant mobility and normalized physical 
quantities) of the Cahn-Hilliard-de Gennes model [9] of phase separation in binary, incompressible, 
isothermal mixtures of polymer molecules of different types i = 1 and i = 2, quenched below a critical 
temperature. 

Following the Cahn-Hilliard-de Gennes theory [5], [TU] (see also [3J, [S], [13), we describe 
now briefly the physical basis of the presented model. Each molecule i in the mixture consists of TV; 
segments of size Oi (so-called lattice constant) with the quantity R 2 . = \Ni(j 2 denoting the mean 
radius of gyration of the i-th polymer molecule. The variables u and w in Problem (P) have the 
meaning of the rescaled order parameter and the exchange chemical potential. For clarity of physical 
description let us use just now in place of u € [—1,1] the original variable X ~ Xi £ [0, 1], denoting the 
volume fraction of component i — 1. Then, by the incompressibility condition, X\ + X2 = 1 everywhere 
in the sample. The variable 

u:=2X-le[-l,l] (1.8) 

is introduced for the sake of mathematical convenience. 

The model is governed by the Flory-Huggins-de Gennes free energy functional [5] which in the 
isothermal case is 

FFHdoiX) = J (F FH {X) + i a (x)|VX| 2 ). (1.9) 
The homogeneous (volumetric) free energy Ffh(X) has the Flory-Huggins form 

Ffh(X) = -^XlogX+-Ul-X)log(l-X) + AX(l-X), (1.10) 

iVi I\ 2 

where A, called the Flory-Huggins interaction parameter, measures the strength of interaction between 
two kinds of species. The second term in the integrand (jl.9l) gives the weighted contribution due to 
composition gradients. By de Gennes theory [S], it has the characteristic singular form 



2 ' ' GXNtX N 2 (l-X)J 36 VX 1-X 

which is due to the connectivity of the chains that constitute the polymer molecules. In the symmet- 
rical case Ni — N 2 — N, o\ = a 2 = 0, R 2 gi = R 2 g2 = R 2 g , the expression (jl.lll) simplifies to 

-a(X) = — — r . (1.12) 

2 v ; 36X(1-X) v ' 

Here we note that the singular form of a(X) introduces an infinite energy penalty near the pure phases. 
This fact turns out to have an important mathematical consequence related to the separation property 
of the solution (see (|2.23p and Remark (|2.6[) below). 

Under the incompressibility condition, X\ + X 2 = 1, the conservation of mass for component 
1 is given by 

X t +V-j=0, j = -A(X)Vw. (1.13) 

Here j is the mass flux defined as the product of the Onsager diffusion coefficient A(X) (effective 
mobility) and the gradient of the exchange potential w which is the difference between the chemical 
potentials of the components, w = w\ — w 2 . Shortly, we shall refer to A and w as the mobility and 
the chemical potential, respectively. 
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According to the derivations in [5], [3], [5], the expression for A takes the form 

A(X) = 7 ^- = X(1-X)A , (1.14) 

where A, = X^Ao is the Onsager coefficients for the i-th component and Ao is a positive constant. 

We point out that the concentration dependence of the mobility A is also typical for phase 
separation of small molecule systems, described by the classical Cahn-Hilliard equation [5], [7]; see 
also [T3], [H], [2D] and the references therein. 

The chemical potential w is defined as the first variation of the functional (|1.9p yielding the 
following equivalent expressions: 



w 



SJ-FHdG 



8X 
fh( 



F,'.„ l \ ) + \'{X)\WX\ 2 - div(a(X)VX) 



F FH (X)-±a'(X)\VX\ 2 -a(X)AX (1.15) 



- F' FH (X) - V^(X)div(^(X)VX). 

As in the standard Cahn-Hilliard theory, the chemical potential can include an additional term eXj, 
with a positive constant e, accounting for possible viscous effects. 

The equations (|1.13p . (j 1 . 14|) . and (jl.151) with the term eXt, e > 0, lead to the following 
degenerate singular Cahn-Hilliard-de Gennes polymer system: 



X t - A div(X(l - X)Vw) = 0, (1.16) 
w = -a(X)AX - ^|VX| 2 + F' FH {X) + eX t (1.17) 



to be considered together with the initial and appropriate boundart conditions. 

As a first step in the study of such system we consider the case of constant mobility by 
linearizing the principal part X(l — X) in (|1.16D . Moreover, for mathematical convenience, we replace 
X by the variable u defined by (|1.8|) , restrict ourselves to the symmetrical mixture (11.121) , and set all 
physical constants equal to unity. This leads to fj 1 . 1 [) - (11 .2[) . 

Apart from the well-known Cahn-Hilliard-de Gennes model described above, we mention also 
recently developed two-fluids models for viscoelastic phase separation in polymer solutions, see [2"9"] 
and the references therein. 

It is known that the process of polymer phase separation is important both for its theoretical 
aspect and due to unusual morphology for specific materials applications. On the contrary to exten- 
sive physical and numerical literature on polymer mixtures (for review see e.g., [3], [ID], [21], [13], [I]), 
mathematical aspects of the corresponding models have not been so far addressed sufficiently. The 
structure of steady-state solutions to degenerate singular polymer system (11.16p - (|1.17[) has been anal- 
ysed by Mitlin et al. [IB] , [H] , and Witelski [35] . The existence of weak solutions to the Cahn-Hilliard 
equation with logarithmic potential F(u), degenerate mobility A(u) = 1 — u 2 , and constant gradient 
coefficient a > has been studied by Elliott and Garcke [12]; we refer also to [5J and [14] for further 
study of degenerate problems. For the standard Cahn-Hilliard equation with logarithmic potential 
F(u) and constant coefficient a > 0, optimal regularity of weak solutions is analyzed in [19) . where in 
particular the separation property (|2.23[) is obtained in space dimensions 1 and 2. The existence and 
uniqueness of weak solutions to the Cahn-Hilliard equation with logarithmic potential and nonlinear 
positive, bounded coefficient a(u) has been recently proved by the authors in [26]. 

We mention also a closely related sixth order Cahn-Hilliard type problem with nonlinear 
coefficient a(u), considered recently in [25] for a singular (e.g., logarithmic) potential, and in [23] for 
a polynomial potential. As a special case, in |26j the behavior of the solutions when the sixth order 
term is let tend to zero was analysed. 

Finally, we point out that the theoretical investigation of polymer models was initiated by Alt 
and the second author in [T], where general nonisothermal phase transition models with a conserved 
order parameter have been derived and, in particular, polymer free energy models have been presented 
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along with an extensive list of references. The authors of [I] have obtained also some partial, unpub- 
lished results [2] on the existence of weak solutions to degenerate singular polymer model (|1.16[) - (|1.17p 
by applying the methods due to Elliott and Garcke [12], and Elliott and Luckhaus [13] . This unsolved 
problem has become the motivation of the present study which uses a different approach developed 
previously in |26l . 

As already mentioned, as a first step of the analysis we assume that the mobility is constant, 
and extend our methods applied in |26] in the case of a nonlinear (but bounded) coefficient a{u) to 
singular a(u). 

There are two main ideas behind our approach. The first one, standard in the analysis of 
Cahn-Hilliard systems, exploits the characteristic variational structure of the system (|l.ll) - (ll.2p . The 
second non-standard one consists in introducing appropriate changes of variables. 

The variational structure becomes evident by (formally) testing by w, (|1.2[) by u t , taking 
the difference of the obtained relations, integrating with respect to space variables, using the no-Hux 
conditions (jl.4[) . and performing suitable integrations by parts. Then one readily gets the a-priori 
bound 

±8(u) + ||VHIi a(n) + e||ut|li*(n) = °- ( L18 ) 
which has the form of an energy equality for the energy functional 

£(u) = |j^M|Vu| 2 + FH-A^), (1.19) 

where the interface (gradient) part contains the nonlinear function a. In other words, the system 
(|1.1[) - (|1.2[) arises as the (7J 1 )'-gradient flow problem for the functional £. 

However, the energy estimate (| 1 . 18[) is not sufficient to obtain existence of a solution to (jl.ip - 
(|1.2j) via approximation-compactness methods. Actually, to apply this strategy, one also needs some 
control of the second space derivatives of u and of the singular coefficients a and a' in (|1.2j) . To obtain 
this, two changes of variables will play an important role. The first one, motivated by the structure 
of the fourth formula in (|1.15[) and applied also in [55], consists in introducing the variable z such 
that Vz = ^y'a^JVu. Then, the Laplacean of z can be (formally) estimated simply by testing the 
equivalent z-formulation of (|1.2[) (namely, (I2.13[) below) by — Az. This gives the desired control on 
second space derivatives. 

On the other hand, even in the equivalent formulation (12 . 13[) . one has to control a coefficient 
(namely, <j/(u)), that explodes polynomially fast as \u\ approaches 1. This is a nontrivial issue since 
the nonlinear term f(u) in (|1.2I) (or in (I2.13[) ), for which is relatively simple to get a L 2 -control, 
explodes only logarithmically fast (and, hence, an L p -estimate of f(u) would help to control (f)'(u), 
or a(u), only for p = oo). To overcome this further difficulty, a second additional change of variable 
comes into help. Namely, we set v := /(«), which represents the monotone part of the volumetric 
chemical potential. We then see that the formulation of equation (ll.2[) in terms of v (namely, (|2.17[) 
below) does no longer contain singular coefficients of polynomial type; it presents, however, the cubic 
term u|Vu| 2 . To control it, we use techniques based on entropy-type estimates (cf. [8]), with a rather 
careful and ad-hoc choice of test functions (cf. (|4.40|) and (I4.63[) below). In this way we can both 
prove an L°°-bound for v (and consequently the "separation property" (|2.23[1 ). for strictly positive 
times, and also control the cubic term it|Vv| 2 starting from the initial time t — 0. This is the key step 
that permits to get existence of a weak solution (for the w-formulation of the system) for initial data 
Mo that have only the natural energy regularity (i.e., such that £ (uq) < +oo). 

Regarding additional properties of solutions, we remark that, as already noticed in [IT], the 
energy (|1.19[) with the coefficient a given by (|1.7[) is convex with respect to u (up to the lower order A- 
perturbation) . This basic property permits to prove parabolic time-regularization properties of weak 
solutions, as well as uniqueness, in a relatively standard way. 

The plan of the paper is as follows. In the next Section[2]we will present the main assumptions 
and give the statement of our main results. In Section [3] we will prove local existence and uniqueness 
of a strong (i.e., lying in a suitable Holder class) solution. The main a-priori estimates needed in the 
proof of global existence will be detailed in the subsequent Section [4] On the basis of these estimates, 
in the final Section [5] we shall show global existence, uniqueness, and time-regularization properties 
of weak solutions. 
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2 Notation and main results 

Let f2 be a smooth bounded domain of R 3 of boundary r, T > a given final time, and let Q :— 
(0, T) x f2. Let H := L 2 (Q), endowed with the standard scalar product (•, •) and the norm || • ||. Let 
also V := if^fi). We identify H with H' so that the chain of continuous embeddings V C H C V 
holds. We indicate by (•, •) the duality between V' and V and by || • \\x the norm in the generic Banach 
space X. We note as A the weak Laplace operator with no-flux boundary conditions, namely 

A:V->V, (Av,z):= [ Vv-Vz, Vv,zeV. (2.1) 

We also set 

W := {z e V : Az e H} = {z£ H 2 (n) : d n z = on T}, (2.2) 

which is a closed subspace of H 2 (SY). In all what follows we shall assume that F, f and a are given, 
respectively, by (1 1 .5f) . (jl.6p and (jl.7p . Moreover, we will assume that uq is an initial datum having 
finite energy £ (cf. (jTTT9j) ) , namely 

E := £(u ) < +oo. (2.3) 
It is worth noting that, since F has the expression (| 1 . 5(1 . the above is equivalent to asking 

- 1 < u < 1 a.e. in 0, a 1/2 (w )Vii € H. (2.4) 

Moreover, due to (|1.7p . if uo satisfies (|2.4p then uq E V, in particular. Letting, for a generic summable 
function v : f2 — > K, denote its spatial mean value, we will also assume that 

m := («o)n = p. J u Q e (-1, 1). (2.5) 

In other words, we cannot admit the case when uo coincides with +1 (or with —1) almost everywhere 
in Q. This is a standard assumption when dealing with Cahn-Hilliard systems containing constraint 
terms (cf, e.g., jTB] for more details). 

Let us note also that, substituting the expression (| 1 . 7[) for a, (|1.2p can be rewritten as 

2 2u 

w = — ^ Aw - ^|Vm| 2 + f(u) -Xu + eu t . (2.6) 

l — u 1 (l — u i y 

It is now convenient to introduce a couple of additional variables which permit to give alternative 
formulations of (|2.6p . To start with, we compute some derivatives of a. From (|1.5p and fj 1 . T[) . we have 



«'W =/"W = (1^2)2 . ( 2 -7) 

as well as 

"M f'»M 4 ( 1 + 3r2 ) ^9 8^ 

a (r) =/ (r) = (1 _ r2)3 ■ (2.8) 

In the sequel, for a generic locally integrable real- valued function tp defined in an open neighbourhood 
of 0, we will write 



tpir) := / ip{s) ds. (2.9) 
Jo 

Then, of course, we have 

a(r)= [ a(s) ds = /(r) = log(l + r) - log(l - r). (2.10) 
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Next, we introduce 

ttr) :=^J Q r a 1/2 ( S ) ds = jT _ d * = arcsinr. (2.11) 
Then, we notice that, setting 

z := 4>{u) — arcsinw, (2-12) 
equation (11.21) can be rewritten in the equivalent form 

w = -2(/)'(u)Az + f(u)-\u + su t . (2.13) 

Next, we put v := f{u). Then, a simple computation gives 

u = r 1 (v) = ^ T \=--j(v). (2.14) 

Moreover, (|1.2j) can be rewritten as 

w = -Av + v+^y^-\Vu\ 2 - Xu + eu t . (2.15) 



Noting that 

" "'V^l 2 = t-^^IV^I 2 = ^\f{u)Vu\\ (2.16) 



2 |V 1 (1-^)2 



we finally obtain from (|2.14l) 

j(v) 1 e° — 1 

w = -Av + v+ ^-^|Vv| 2 - Xj(v) + eu t = -Av + v H |Vv| 2 - Xj(v) + eu t . (2.17) 

2 2 e v + 1 

In the sequel, we shall indicate by c a generic positive constant, whose value may vary on 
occurrence, allowed to depend on the parameters of the system (more precisely, on the functions a 
and /, on A and on f2), and, in particular, not on approximating parameters. Moreover, the constants 
c will not be allowed to depend on the choice of initial data. However, they may depend on the 
prescribed mean value m. The notation k will be used for positive constants (depending on the same 
quantities as c) appearing in estimates from below. We will also use the notation Q(-) (or <5(-,-)), 
with Q indicating a computable function with values in [0, +oo), increasingly monotone in each of 
its argument, whose expression can depend on the same quantities as c. For instance, the expression 
Q(Eo) will stand for a monotone function of the "initial energy" Eo- 

We can now introduce the concepts of "classical" and of "weak" solution needed in the sub- 
sequent analysis. 

Definition 2.1. A "strong", or "classical", solution to Problem (P) over the time interval (0,T) is a 
couple (u, w) with the regularity 

ue VK 1 ^ oo (0,T;T/')ni? 1 (0,T;t/)n J L oo (0,T;iJ 2 (^)), eu € W /1,oo (0, T; H), (2.18) 
w e L°°(0,T;V)nL 2 {0,T;H 3 (n)), ew G L°°(0, T; H 2 {Vt)), (2.19) 

satisfying, a.e. in (0,T), the equations 

u f -Ara = 0, a.e. in CI, (2.20) 

w = -a(u)Au- ^-t^|Vu| 2 + f(u) - Xu + eut, a.e. in fi, (2.21) 
d n u = d n w = 0, a.e. on T, (2.22) 

together with the initial condition u\t=o = uo and, for all (t,x) G [0, T] x fl, the separation property 

- 1 + e < u(t, x) < 1 - e, for some e > 0. (2.23) 



G 



Remark 2.2. Thanks to (|2.23p . the component u of any "classical" solution to Problem (P) is 
uniformly separated from the singular values ±1 of a and /. Hence, by applying the standard theory 
of quasilinear parabolic equation and a bootstrap argument, we can see that (u, uu) is in fact smoother, 
with its regularity being limited only by the regularity of the initial datum. In other words, at least 
for times t > 0, a classical solution can be thought to be arbitrarily regular. 

Definition 2.3. A "weak", or "energy", solution to Problem (P) over the time interval (0,T) is a 



couple (u, w) with the regularity 

u G H x (0,T; V) nL°°(0,T; V) nZ°°((0,T) x Q), eu G H 1 (0, T; H), (2.24) 

F{u) &L°°(Q,T;L l {Ci)), (2.25) 

v = f(u) e L 2 (0,T;V), Vu G L p ((0,T) x Q) for some p > 2, (2.26) 

weL 2 (0,T;V), (2.27) 

satisfying, a.e. in (0, T), the equations 

u t +Aw = 0, in V, (2.28) 

w = -Av + v + ^|Vu| 2 - Xu + eu t , a.e. in Q, (2.29) 

v = f{u), a.e. in Q, (2.30) 

d n v = 0, a.e. on T. (2.31) 

together with the initial condition u\t=o = Uq. 



We can now state our main result, regarding existence, uniqueness, and regularization prop- 
erties of weak solutions to Problem (P): 

Theorem 2.4. Let f and a be given by (fTB)) , ([TTT]) . Let u satisfy (1231) and (|23|) . Finally, let il be 
convex. Then, for any T > 0, Problem (P) admits at least a weak solution (u, w) dehned over (0, T). 
Moreover, for any r > 0, (it, w) is a "classical" solution over (r, T). In particular, the separation 
property ([2^3]) holds on (r, T) with 

e- 1 = Q(E ,t- 1 ). (2.32) 
Finally, uniqueness holds in the class of weak solutions that are classical for strictly positive times. 

For strictly positive times, (|1.2p can be interpreted in any of the equivalent formulations (|2.6[) . (|2.13l) . 
or (|2.17p : actually, u is a classical solution for t > 0. On the other hand, when looking at the behavior 
near t = 0, it is crucial to view (11.21) in the form (12 .29[) which appears to be the only formulation 
permitting to take the approximation limit S \ starting from the initial time. 

Remark 2.5. If we have in addition that the initial datum satisfies uq £ H 2 (fl) with uq(x) G 
[— 1 + e.Q, 1 — eo] for all x G fi and some eo G (0, 1), there are no complications due to the boundary 
layer t \ 0, and so u can be seen as a "classical" solution over the whole (0, T). This can be deduced 
simply by using the classical Gronwall lemma (instead of the uniform Gronwall lemma) in the a-priori 
estimates. As a consequence, estimates (|4.30p . (14.601) and (|4.62l) below hold in fact with t = in this 
case. 

Remark 2.6. It is worth stressing that the separation property (|2.23l) holding for our model is instead 
an open issue in the frame of the standard Cahn-Hilliard system with logarithmic nonlinearity f(u), 
at least in the three-dimensional setting (cf. [TH] for further remarks). 

3 Local strong solutions 

In this section, we will prove existence of at least one local in time classical solution to Problem (P). 
With this aim, we first introduce a regularization of the initial datum uq. This is the object of the 
following 
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Lemma 3.1. Let uq satisfy (|2.3[) and (|2.5p . Then, for any 8 G (0, 1/6), there exists uq } s such that 
u Q j G C°' a (fi) for any a G (0, 1/2). Moreover, 

- 1 + 36 < u 0t s{x) < 1- 3<5 VxGfL (3.1) 

Finally, we have that 

uo,S uo weakly in V, (3-2) 

£(u ,g) <c(l+£(u )) V<5g (0,1/6). (3.3) 



Proof. First of all, we set 



1 o.s 



min {l — 3(5, max{u , — 1 + 3(5}}, (3.4) 



so that Mq 1 ] satisfies a.e. in Q the equivalent of Q3.ip . Moreover, it is clear that £ (itg 1 ]) < £{uq) for all 

6. Next, we define Zg 1 ] := (^(u^l) — arcsiniig 1 ]. Then, we proceed by singular perturbation, defining 
zg ,5 as the unique solution of the elliptic problem 

zo,s + 8Azo,s = Zqj- ( 3 - 5 ) 

Being Zq 1 ] G V, then, by elliptic regularity (recall that Q is a smooth domain), zqj G i/ 3 (f2) for all 6. 
Setting uo,g ■— sinzo.a, a direct check permits to verify that (at least) uqj G H 2 (Q) for all S. Hence, 
ug,5 is Holder continuous, as desired. Moreover, by monotonicity of <f> and a standard maximum 
principle argument, it is clear that (|3.1|i holds. 

The key step consists in proving (|3.3[) . Actually, it is obvious that ||-F , (t«o,5)||ii(n) < c 
(cf. (|2.4p V To control the gradient term of £ , we test (|3. 5[) by Azq,<5- We obtain 



aK«5)|V Uo> | 2 = ||Vz 0i(5 || 2 < ||VzJj|| 2 



as desired. Hence, we have (13.3p . To conclude, we have to prove (|3.2p . which is however an immediate 
consequence of standard weak compactness arguments. The proof is complete. I 



a(u$)|Vu$| J < - / a(u )\Vu \ 2 < £(u ) + c, (3.6) 



As a next step, we also provide a modification of the function / given by (|1.6p . Namely, for all 

6 G (0, 1/6), we take fg G C l {-1 + 6, 1 - 8) -> R such that f$ is monotone, |/ 5 (r)| > for all 
r G (-1 + 5, 1 - (5), and 

fs(r)=f(r) if re [-1 + 2,5,1-2(5], (3.7) 

lim fs(r) signr = +oo. (3-8) 

|r|— 6 

It is obvious that, for any <5 G (0, 1/6), a function fg with the above properties exists. 
Finally, we modify a by taking ag G C 2 (R;K) such that 

a s (r)=a{r) Vr G [—1 +8,1 — 6]. (3.9) 

In particular, 

ojf(r)>0, (— )"(r) = -l VrG [-1 + 8,1-8]. (3.10) 
\ag/ 

Outside (—1, 1), as is taken as a constant Kg > (exploding as 8 \ 0), whereas for |r| G (1 — 8, 1), 
ag is chosen in such a way to have 

l<ag(r)<Kg VreR. (3.11) 
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Then, for 5 £ (0, 1/6), we can consider the system 



u s ,t-Aws=0, in(0,T)xQ, (3.12) 

w s = -a s (u s )Aus - ^y^|Vu A -| 2 + f s (u 5 ) - \u 6 + eu s . u in (0,T) x Q, (3.13) 

ug\ t =o = u j, in ft, (3-14) 

d n us = d n ws — 0, on T. (3.15) 

We then have: 



Theorem 3.2. Let f and a be given respectively by (jl.6l) . (11.71) . and let uq satisfy p.3[) and (|2.5I) . 
For <5 £ (0, 1/6), Jet uo,i5 be defined by Lemma |3~T1 and fs, as be given by (|3.7p - (|3.11[) . Then, there 
exists one and only one solution (us,ws) to system (|3.12[) - (|3.15[) with the regularity 

u s £ W 1 ' 00 ^, T; 7') n H^O, T; 7) n r (0, T; 2 (fi)), eu s £ W h °°(0, T; H), (3.16) 
w 5 £L°°(0,T;V)nL 2 (0,T;H 3 (n)), ew s £ L°°(0, T; H 2 (Q)). (3.17) 

Proof. We claim that this result is essentially a consequence of the results of [2BJ. Actually, we 
see that, for any 5 £ (0,1/6), as satisfies the assumptions [25] (2.1)-(2.2)] (where —1 and 1 replace 
-2 and 2 in [S3 (2.2)]) and US (6.1)] (where [-1, 1] is replaced by [-1 + 6,1- 5], cf. (l3~TT)l) 'l. 

Moreover, satisfies [26j (2.3)-(2.4)], with (—1 + 5, 1 — S) replacing (—1, 1). Thus, we can 
apply 26, Theorems 5.1, 6.1, 6.2] which give the existence and uniqueness of a weak solution (us,wg) 
to (|332>(|3~1"5]) . The regularity of this solution is specified by [23 (4.3)-(4.4) and (6.14)]. More 
precisely, since the initial datum uq.$ is smooth and separated in the uniform norm from the singular 
values ±(1 — S) of fs due to (|3.1[) . we have here that [23 (6.14)] holds starting from the initial 
time, i.e., with r = 0. Moreover, a closer inspection of [251 Proof of Theorem 6.2] (see in particular 
estimates (6.19)-(6.20) therein) permits to see that the additional regularity for Us,t stated in (|3.16[) 
holds over (0, T). The L 2 (0, T; H 3 (n)) regularity of w s follows from the L 2 (0, T; V) regularity of u S j 
and elliptic regularity estimates applied to (|3.12p . 

Hence, collecting all the information coming from the results of [25], we obtain exactly (|3.16D - 
(|3.17[) . This concludes the proof. I 

Notice now that, as a consequence of (I3.16[) - (|3.17[) and of the arguments in [55], we have, more 
precisely, the a-priori estimate 

li^aliffHo, t-v) + \\ u s\\L°°{a,T;H 2 (n)) < Q{\\uo,s\\h 2 (q)>^ )• (3.18) 
By interpolation and embedding properties of Sobolev spaces, we then obtain 

H M <5llc < b ([o,T]xO) ^ QiWuoMlH^n),^ 1 ), forsome6>0. (3.19) 

Thus, us is globally Holder continuous. Since iio,a satisfies (|3.1[) . we can deduce that, once the initial 
datum uq is assigned, then for all 5 £ (0, 1/6) there exists a (computable) time Ts > depending on 
S and uq such that 

- 1 + 25 < u s (t, x)<l-26 for all (t, x) £ [0, T 5 ] x Ti. (3.20) 

However, over [— 1 + 26, 1 — 2(5], as coincides with a by (|3.9[) and fs coincides with / by (13.71) . Hence, 
we have obtained the 

Corollary 3.3. Let / and a be given by CE1|, $1J§, and let u satisfy (23]) and ([23]). Let 5 £ (0, 1/6) 
and let Uo t s be given by Lemma |3. II Then, there exist a time Ts > depending on uq and S, and a 
"classical" solution (us,u>s) to Problem (P), with initial datum uqj, over the time interval (0,Ts). 



9 



4 A priori estimates 

In this section, we derive a number of a-priori estimates for the solutions of our system. We stress that 
the procedure leading to these estimates can be rigorously justified at least for "classical" solutions 
to Problem (P) in the sense of Definition 12.11 Indeed, owing to Remark 12.21 these solutions can be 
thought to be as smooth as wc need (possibly paying the price of additionally regularizing uq). 

In particular, the estimates proved below will hold for local strong solutions given by The- 
orem 13.21 On the other hand, for weaker notions of solutions the procedure below may just have a 
formal character due to insufficient regularity of test functions. We will clarify this point in Section [5] 
below. Here, we will proceed assuming that everything is regular enough for our purposes. Moreover, 
since the estimates we are going to derive will have a global-in time character, with some abuse of 
notation we will directly work on the time interval [0,T]. The underlying extension argument will be 
also detailed in the next section. 

We can now start detailing our estimates. 

Energy estimate. We test (jl.ll) by w and (| 1 . 2|) by u t . This gives rise to the energy equality (jl.181) . 
whence we obtain the estimate 

IMIl2(0,T;V) + £ \\ U t\\h(0,T;H) + II ^ W \\h{Q.T;H) < Q(^o), (4-1) 

II^WIU~(0,T;iMO)) < COM- (4.2) 
Thus, by definition (|1.5[) of F, we infer in particular that 

-1<«<1 a.e. in (0, T) x fi. (4.3) 
Moreover, recalling also (12.12[) . we obtain 

ll«IU»(o,r 5 v) + ll z IU~(o,T;V) < <9(Eo). (4.4) 

Estimate on time derivatives. Following the lines of [551 Theorem 6.1], we indicate by J the 
gradient part of the energy, i.e., 

J:V->[0,+oo), J{u) := f ^\Wu\ 2 . (4.5) 

Then, we can (formally) compute the first derivative of J, given by 

a'{u) 



(J'(u),v) = j (a(u)Vu-Vv+?-^\Vu\ 2 vy (4.6) 
as well as the second derivative 

(J"(u)v, z)= j ( a "( u )\J u \ 2vz + a '{u) v ^u . Vz + a '(u)zVu ■ Vv + a(u)Vv ■ Vz) . (4.7) 



To be more precise, if u is a smooth solution (and in particular it is separated from singularities, i.e., 
it satisfies (|2.23l) ). formulas (|4.6[) and (|4.7I) make sense for u = u(t) at any time t £ [0, T] and, indeed, 
wc have that J'{u) € V and J"{u) € C(V, V). 

From (|4.7p . we then obtain in particular 



J"(u)v,v)= / ^ 2 + 2a'{u)vVu- V v + a{u)\V v\ 2 j 



n -> 



> / W~ 2 4£t-)\Vv\ 2 . (4.8) 



n x a"(u) , 

We can now test (|1.1[) by twt and add the time derivative of (|1.2p tested by tut- This leads to 

l^ 11 VH|2 + ld7 l|Uil ' 2 + 1 (J>K Ut) + 1 J n f,(u)u2t = At||utl ' 2 - (4 - 9) 
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Hence, using (|4.8[) . noting that 

2a'(u) 2 



°^-W = TW (4 - 10) 



and using (|4.3I) . we arrive at the inequality 



^^IIVHI 2 + f INI 2 ) + 1 J + */ !'iu)u 



< Atll^ll 2 + i||V W || 2 + £ -\\u t \\ 2 < t-WVutf + c(l + t)||V^|| 2 + ||K|| 2 . (4.11) 

Recalling (|4~T|) and noting that 2/(1 + 3-u 2 ) > 1/2 since u takes values in [—1, 1], integration in time 
of (|4"TT|) gives, for any t <E (0,T), 

+ || Vlfl||£- (T ,T;JO < Q( E 0> ( 4 - 12 ) 

e||«*||l-(t,t ; h) + H M *IU 2 (r,T ; y) < Q^O. (4.13) 
Estimate of /(it). We test (|1.2[) by u — m. Integrating by parts the terms depending on a, we obtain 

(u) + a ^ (u - m)\ \Vu\ 2 + / f(u)(u — m) = (w + Xu — su t ,u — m). (4.14) 

We have to estimate some terms. Firstly, proceeding as in [TIB Appendix], it is not difficult to prove 
that _ 

- m) > -\\f(u)\\ L i (n) - c m . (4.15) 

We notice that assumption (I2.5[) is used here. 
We now observe that 

/„ («(») + - ■»>) i'»i 2 ^ /„ + rh) <««> 

Moreover, there exist constants K m > 0, c m > such that 

2u(u — m) 1 u 2 . . 

(1 - U 2 ) 2 ^ K "» (1 _ U 2)2 - C ™ ^ K ™ (1 _ M 2)2 - C «- ( 4 - 17 ) 

Thus, recalling (|1.7[) . we get 

(o(«) + ^-{u- m)) \Vu\ 2 > K m \\Vf(u)\\ 2 ~ c m \\Vu\\ 2 . (4.18) 

Finally, noting that u takes values in [—1,1], using estimate (I4.4[) . and applying the Poincare-Wirtinger 
inequality, we have 

(w + Xu — eut, u — m) — (w — tuji + Xu — eit t , u — to) < c(|| Vw| + 1 + e||ut||) ■ (4-19) 
Thus, collecting the above considerations, (|4.14[) gives 

||V/( U )|| 2 + ||/(«)|Ui ( n) < c(ll V «l| 2 + H V HI + 1 + e IM)- (4.20) 
Squaring (|4.20p . using <|4. 15|) - (|4. 19[) . and integrating in time, we arrive at 

l|V/( U )|| L 4 (0iT;H ) + ||/(u)||i 3 (o,Ti£i(n)) < Q(Eo), (4-21) 

whence, more precisely, 

\\f(u)\\mo,T;V) < Q(Eo). (4-22) 

Taking instead the essential supremum of (|4.20p as t ranges in (r, T), and using (|4.12p - ([4.13l) . a 
straighforward modification of the above procedure leads to 

ll/(«)IU»(r,T;V) ^Q^o.r- 1 ) Vre(0,T). (4.23) 
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Next, integrating (jl.2l) in space, using (|2.7[) . (|4.3|) . (|1.7|l . and noting that (ut)n = 0, we get 

\w n \ < c( jf ^|Vu| 2 + ||/(«)|U 1(n) + l) 

< C (||V/( U )|| 2 + ||/( U )|| il(0) + l). (4.24) 
Thus, squaring, integrating in time, using (|4.2ip . and recalling the last (|4.ip . we infer 

H\lho,T;V) < 0(Eo). (4-25) 
Taking the essential supremum in (|4.24D as i <E (r, T) , and recalling (|4.12[) and (|4.23p , we also get 

IMU-frTiV) < 0(^0,0 Vre(0,T). (4.26) 
Estimate of z. We consider the equivalent formulation (|2.13p and test it by — Az. This gives 

/ 2(j)'{u)\Az\ 2 + (f'(u)Vu,Vz) = (w + Xu-eu t ,-Az). (4.27) 
Jn 

Thus, using the monotonicity of / and <p, noting that 4>'(u) > 1/2 for all u G (—1,1), and that 
(/'(u)Vu, Vz) > 2|| Vz|| 2 , we can control the right hand side this way: 

(w + Xu-eu u -Az) < i||Az|| 2 + c(||w|| 2 + l + e 2 ||7i t || 2 ). (4.28) 
Hence, integrating (I4.27[) in time and recalling (|4.ip , (|4.4I) and (I4.25[) . we arrive at 

||«||i»(o,Tjfla(n)) < Q(Eo). (4-29) 

Taking instead the (essential) supremum of (|4.27p as t ranges in (r, T) for r > and using (|4.12l) - (|4.13p 
and (|4.26l) . we obtain 

IMU~(r,T;fl»(n)) < Q^o,^ 1 ). (4.30) 

The above relations permit to improve also the bounds on u. Actually, computing directly the 
Laplacean of u = sinz and using f)4.3[) together with the Gagliardo-Nirenberg inequality (cf., e.g., 
[2"2"1 Theorem p. 125]) 

l|Vy|U* (n) < co||y||^ 2 (0) ||2/||^ (0) + ||y|| Vy e H 2 (V), (4.31) 
it is not difficult to arrive at 

\\u\\mo,T-m(n)) <<3(Eo), (4.32) 
||«IU-(r,T;fl»(n)) < QQEo,!-- 1 ). (4.33) 

First entropy estimate and separation property. The estimates obtained up to this moment 
yield a control of the functions z and u up to their second space derivatives (cf. (|4.29[) - (14.301) and 
(|4.32l) - (|4.33l0 . and of the nonlinear term f{u) up to its first space derivatives (cf. (|4.22[) and (|4.23t ). 
However, this still seems not sufficient to pass to the limit in the equation (|1.2I) . even if its equivalent 
formulation (|2.13p is considered. Indeed, from (|4.22p and (I4.23P we get a control of the term f(u), that 
explodes logarithmically fast as |u| /*• 1. On the other hand, even in formulation (|2 . 1 3[) one faces the 
term (j)'(u) which is much more singular since it explodes as a negative power of (1 — u 2 ). To control 
it in some L p -norm we need more refined estimates of the so-called entropy type and, in particular, 
we need to refer to the formulation (|2 . 1 7[) (we recall that all formulations are equivalent, at least for 
sufficiently smooth solutions). Usage of this technique requires the convexity assumption on f2 asked 
in the statement of Theorem 12.41 

The basic tool we need consists in an integration by parts formula due to Dal Passo, Garcke 
and Grim ((5J Lemma 2.3]): 
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Lemma 4.1. Let h e W 2 >°°(R) and zeW. Then, 

f h'(z)\Vz\ 2 Az = -\ f h"(z)\Vz\* 
Jn 6 Jn 

+ \J h(z)(\D 2 z\ 2 - \Az\ 2 ) + | J h(z)II(Vz), (4.34) 

where II(-) denotes the second fundamental form ofT. 
Then, we test (j!TL7f by 

— Am(w) = — div(m(v)Vv) = —m(v)Av — m' (v)\Vv\ 2 , (4.35) 
with the function m to be chosen later. This gives 

j m[v)\W\ 2 + J (rn{v)\Av\ 2 - ^M|V,| 4 + (m» - ^MM) Av] y v A 

= (m(w)Vf,Vw + AVu) + (eu t ,m(v)Av + m! (v)\Vv\ 2 ) . (4.36) 

Now, applying Lemma 14.11 to the last integral on the left hand side of (|4.36[) we infer 

m{v)j(v) \ 2 Iff „ m'(v)j(v) m(v)f(v) \ 4 
m(v) \Av\Wv\ =--/ m (v) )|Vu| 



2 ) 1 1 3 J n V ' ' 2 2 

3 



- ™ + K) (\D 2 v\ 2 - \Av\ 2 ) + | / (m(v) - + tf) JJ(V«), (4.37) 



where K > is an integration constant that will be chosen later on, and the notation (|2.9[) is used. 

Substituting (|4.37l) into (|4.36[) . we get on the left hand side the following "hopefully good" 

terms 

1 f ( ^nr..\ ^'r„.\Af„.\ , rn{v)j'(vy 



where we notice that 



m(v)\Wv\ 2 + / m{v)\Av\ 2 + - / ( - m"(v) - m'(v)j(v) + mWj W ) |V^| 4 , (4.38) 

.In 3 ./o v 2 / 



We can now specify our choice of m as 



j'M = - 2e (4.39) 

J w (e" + 1) 2 y J 



m(v) = ^—ry (4.40) 



Actually, this expression arises since we need m to decay not too fast at infinity (otherwise we do not 
get enough information from it), but at the same time we need it to be summable (cf. (I4.47[) below). 
The above choice gives 

m(v) = — arctanw, (4-41) 

as well as 

-» = -^TW' m " {v) = (TTW- (442) 

Now, noting that mf > 0, we can easily observe that 

/// \ it \ ' t \ m(v)f(v) ^ (v + v 3 )j(v) + 1 - 3v 2 , A ln . 

-m"(v)-m'(v)j(v)+ W -/ > ^ (1+^2)3 ' ( 4 ' 43 ) 

Now, recalling (|2.14p . for a suitable M > (e.g., we can take M = 12 here), a direct computation 
shows that 

H^) + l-3^11 + tf VH ^ M 



(1 + V 2)3 -4(1 + W 2 ) 3 
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On the other hand, 

(v + v 3 )j(v) + l~3v 2 1 3v 2 l + \v\ 3 \v\ UA ,, 

T« 7T5 > T, nTT — T. oTT > KM7Z ^7 — CM TZ ^TT V \v\ < M. (4.45) 

(1 + u 2 ) 3 _ (l + i; 2 ) 3 (l + i; 2 ) 3- (l + i; 2 ) 3 (1 + d 2 ) 3 1 i_ v ; 
Summarizing, we have 

a, s i , \ , \ m(v)f(v) 1 + M 3 \v\ .„ 

m"(v) m'(v)j(v) + » > ^^^p ^ c^-LL_ X{|t)| < M} , (4.46) 

where X denotes the characteristic function. 

Now, we can choose K so large that the function 

( m( ,)_™M + if ) (4 .4 7) 

is strictly positive (and bounded, of course) . Then, thanks the convexity assumption on f2 (that entails 
positive definiteness of the second fundamental form), the latter two terms in (|4.37j) are positive. 
Collecting these observations, we can deduce from (|4.36[) the estimate 

/ m{v)\Vv\ 2 + [ m(v)\Av\ 2 + n f 1±M1| V «| 4 <E T - ( 4 - 48 ) 
Jsi J si Jsi \ L + v ) i=1 

and we have to control the "bad" terms on the right hand side: 

Ti-=cf H 3 |V^| 4 , (4.49) 
J{\v\<M} I 1 + v ) 

T 2 := / m{v)Vv ■ Vio, (4.50) 
Jn 

T 3 := A / m(t))Vi) • Vu, (4.51) 

J SI 

T 4 := (eu t ,m(v)Av + to»|Vv| 2 ). (4.52) 

To do this, we first notice that 

\T 2 +T Z \<- ( m(w)|Vt;| 2 +c(||Vw|| 2 + ||Vw|| 2 ). (4.53) 
4 J si 

Next, we have to control T\. Then, we can note that there exists 8 £ (0, 1), depending only on M, 
such that the restriction / : [— 1 + 8, 1 — 8] —> [—M,M] is bijective and Lipschitz continuous together 
with its inverse j. Thus, using (I4.3[) and inequality (|4.31[) . we deduce 



/» 4 |Vu| 4 

'{\v\<M} J{H<M} 



\Ti\<c[ \Vv\ 4 <c[ 

J{\v\<M} J{\v\ 

<cf h \w \Vu\ 4 <c 5 I |Vu| 4 

J{\u\<l-8} \ l ~ u ) J{\u\<l-S} 

<c s [ |Vu| 4 < c*||«|||« (n) (l + \\uf H2(n) ) < cj(l + \\uf H2{n) ). (4.54) 
Jsi 

Finally, we have to control T4. We have 

|T 4 | < \ f m(v)\Av\ 2 + c a e 2 \\u t \\ 2 + a f m'{v) 2 \Vv\\ (4.55) 
4 Jsi Jsi 

for small a > to be chosen later. Then, we can go on as follows: 

o [ m '(v) 2 \Vv\ 4 = a [ ^^| Vt ;| 4 <af ±±^-\Vv\\ (4.56) 
Jsi Jsi(l + v 2 r Jsi{i + v 2 y 
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Hence, for a sufficiently small, the last integral is controlled by the last term on the left hand side 
of (|4.48p . On the other hand, the first term on the right hand side of (|4.55p is controlled by the 
corresponding one on the left hand side of (|4.48[) . 

Now, a direct computation permits to see that 

I * i Pi fen 't 1/2 l^| 4 > « I I V(l + - 2 ) 1/8 | 4 - (4-57) 
Thus, collecting ([233 > ([237 | . (jOgjl gives 

m(w)|Vw| 2 + / to(v)|Ai;| 2 + / |V(l + v 2 ) 1/8 | 4 
i Ja 

< C (||Vu,|| 2 + ||V U || 2 + 1 + \\uf H2{n) +e 2 \\u t \\ 2 ). (4.58) 
Hence, integrating in time, and using (14.11) . (|4.4[) and (|4.32p . we obtain 

On the other hand, taking the essential supremum in (|4.58p as t ranges in (r, T) for r > 0, and using 
(j4~T2"l) . |gJ3D and (ginSll , we arrive at 

H^A^^^ + |(1 + « 2 ) 1/8 | iO o (TiT;VK1 ,4 (0)) < Q(Eo, r- 1 ). (4.60) 

By the continuous embedding W 1,4 (f2) C L°°(f2) we have in particular 

|(l + « 2 ) 1/8 | ioo((T , T)xf2) <Q(IEo,r- 1 ). (4.61) 

In terms of u the above estimate gives rise to the separation property 

- l + e< u(x,t) < 1-e a.e. in (r, T) x n, (4.62) 

for all t > 0, with e > depending on r. 

Refined entropy estimate. We repeat the entropy estimate of before taking now, in place of (|4.40p . 

m(v) = -j- 1 9N , (4.63) 
v ; 2(1 + w 2 )p 

where the choice of p £ (1/2, 1] will be made precise later on. Then, we have 

// \ P v in \ (2p 2 +p)v 2 -p 

m («) = - {1 + V 2 )P+ 1 ' m W = {1+v 2 )P+ 2 ■ ( 4 - 64 ) 

Thus, a straighforward modification of (|4.43l) - (|4.45l) leads to 

"l \ V.WU^WWs l + \v\ 3 H , . 

- m (v) - m + > « (1 + t;2)p+2 - c (1 + t;2)p+2 *{M<M}, (4.65) 



whereas the equivalent of (|4.57p gives now rise to 



(4.66) 



Now, let us observe that, since p > 1/2, then the function mj is bounded. Thus, we can still take K 
so large, depending of course on p, that the function in (|4.47l) is strictly positive. Moreover, the terms 
corresponding to Tj, j = 1, ... ,4, can be controlled similarly as before. Thus, integrating in time the 
p-analogue of (|4.48|) . we obtain 

T 

V(l + t ; 2 )t-«r<Q(E ). (4.67) 



JQ 
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Let us now test (|2.17p by v. This gives 

l|Vv|| 2 + |M| 2 + / ^^|Vw| 2 < (w + Xu-eu t ,v). (4.68) 



Then, integrating in time, noting that 1 + j(v)v > k(1 + v 2 ) 1 / 2 , recalling (|4.1I) and using Holder's and 
Young's inequalities to estimate the right hand side, we arrive at 

T 

l + « 2 ) 1/2 |Vtf < Q(E„). (4.69) 



o Jn 
Then, let us define 

n+{t) {x eft : \v(x,t)\>l}, n-(t) := {x G H : |u(a;, f)| < l}. (4.70) 

Let us now see that (|4.67p and (|4.69p permit to prove higher intcgrability properties of Vu. Firstly, 
recalling (|4.4I) and ()4.32[) . and using inequality (I4.31[) . we obtain 

IMU^o.zw^cn)) - Q( E o)- (4.71) 

Hence, noting that j is locally Lipschitz continuous with its inverse, we get an analogous information 
for v in the space-time set where it is small: 

f |Vz;| 4 < Q(Eo). (4.72) 
o Jn-(t) 

On the other hand, in f2+(/j) we can write, for 77 > 0, 

iv V | 2 = (v(i + v 2 r r >vv) ■ ( {1+y2) \ v 



v 

= Ct? V(l + « 2 ) 1 -'"- ((l±plvu). (4.73) 

Thus, choosing 77 = | + |, we have 1 — 77 = § — |. Then, we can take p = \ + e, with e > as small 
as we want. We then obtain that 77 — \ = \ + | . Thus, recalling (|4.67p and (|4.69[) . we arrive at 



|Vw| 2 < c v 



7 (1 + v 2 ) 1 ^ • ((1 + v 2 )^ 4 Vv) + - 



(4.74) 



L 4 L 2 

where we can take r strictly greater than 4 since for |u| > 1 we have 

\v\ \v\ 

and we know from (|4.22l) that v is controlled in L 2 (0,T; H). Thus, estimating the right hand side of 
(|4.74p by Young's inequality, integrating first over Q + (t) and then for t € (0, T), we obtain that 

/ / |Vw| 9 < Q(Eo), for some <? strictly larger than 2, (4.76) 

and, of course, combining with (|4.72[) . 

/ / \^v\ q < Q(E ), for some g strictly larger than 2. (4.77) 
Jo Jn 
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5 Existence and uniqueness of weak solutions 

We detail here the proof of Theorem l2.4l which is largely based on the estimates derived in the previous 
section. As a first step, however, we show uniqueness, which works similarly to [26]. Indeed, the key 
assumption (6.1)] is satisfied by our function a (we have, indeed, that (1/a)" = —1). 

Then, let us take a couple of weak solutions (ui,wi) and (1*2,^2) originating from the same 
initial datum uq, and assume that both are "classical" (and in particular satisfy the separation property 
(|2.23[) ). for strictly positive times. 

Setting (u, w) := (ui,Wi) — (u 2 ,W2), we can write both (|2.28[) and (|2.29[) for the two solutions 
and take the difference. Using notation (|4.5p we get 

u t + Aw = 0, (5.1) 
w = J'{ux) - J'{u 2 ) + /(ui) - f(u 2 ) - Au + eu t . (5.2) 

Then, we can test (|5.1[) by A^u, (|5.2[) by u, and take the difference. Actually, the operator A is 
invertible as it is restricted to 0-mean valued functions (as in the case of u due to conservation of 
mass). Noting that 

(Aw, A~ u) = (A(u> - w n ), j4 -1 u) = (w - wu,u) = (w, u), (5.3) 

we then obtain 

5^(lN 2 v" +4u\\ 2 ) + (J'(ui) - J'(u a ),u> + (f( Ul ) - /(«a),u) < A|M| 2 . (5.4) 
Then, recalling (|4.8p and (|4.10[) and using monotonicity of /, we arrive at 

~(\\u\\l r +e\\u\\ 2 )+^\\Vu\\ 2 <\\\u\\ 2 . (5.5) 
Noting that, by the Poincare-Wirtinger inequality, 

A|M| 2 < i||V U || 2 + C |M| 2 y „ (5.6) 
we can integrate (|5.5p over (r, T) for r > 0. Using Gronwall's lemma, we obtain 

||«i - "2||ioo (T)T; ^) < c(T)\\ Ui {t) - u 2 {t)\\ 2 v ,, (5.7) 

where c(T) is independent of r. Then, uniqueness follows by taking the limit r \ and owing to 
continuity of weak solutions with values in V (which is an obvious consequence of (12.2411 ). 

Let us now switch to existence. To start with, we approximate the initial datum uq as specified 
in Lemma \3. II Then, thanks to Corollary 13.31 for any 5 £ (0, 1/6), there exists a "classical" solution 
(us,ws) to Problem (P) defined at least on the time interval (0, Tg), where Tg depends on ug and S. 
Actually, in principle, we may have that T$ \ as we let 6 \ 0. On the other hand, the forthcoming 
argument will exclude this eventuality and show that, in fact, (ug, Wg) can be extended up to the final 
time T . 

Indeed, let us denote as Tg^ max the maximum time up to which (ug,wg) can be extended in 
the form of a "classical" solution; namely, 

T^max := supjS 1 £ (0, T] : ug admits a "classical" extension over (0,5)}. (5.8) 

Due to uniqueness proved above, all extensions of (ug,wg) can be "glued" together. Consequently, 
there exists a (unique) maximal classical extension (ug <max , wg tinax ) defined over (0, T,5 )max ). We claim 
that X^max = T, and, to prove this claim, we proceed as usual by contradiction. Actually, due to 
(|2T8| - (jQ9)) and ([2~^3jl . for any 5 £ (0, T s , ma , x ) we have 

||u < 5,max||Hi(0,S;^)nW 1 -° c (0,S;V)ni° c (0,S;H 2 (O)) <C(5), s\\ M^ max || W 1 '™ (0,S;H) <C(5), (5.9) 

||W5, max |j L oo( 0iS .y) nL 2( 0S . ff 3( n )) < C(5), £||W5, max || L oo( 0:5 . H 2( f2 )) < C(5), (5.10) 

-l + e(S)<u a , ni «(t > i)<l-e(S) V (t, x) £ [0, 5] x Q, (5.11) 
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where C(S),e(S) > and it may be C(S) /> oo and e(S) \ as 5 / T«, m <uc< On the other hand, 
since (ug. max 5 1^8, max 

) is a "classical" solution, it satisfies the a-priori estimates of the previous section 
on the time interval (0,Ts tmax ). Then, thanks to (|4.62l) . we have that 

-l + e<u s . max {t,x) < 1-e V (t, x) £ [t, S\ X ft, (5.12) 

and for all < t < S < Ts, max , with e independent both of S and of S. To be more precise, we have 
that 

e- 1 = Q(£(u , 5 ), t- 1 ) < Q(E , r- 1 ), (5.13) 

where the second inequality is a consequence of fl3.3[) . 

Analogously, we have estimates of the norms in (|5.9[) - ([5.10p over the time interval (r, S) by a 
constant C independent of S and S. Consequently, we obtain 

3u= lim u s , meix (t,-). (5-14) 

To be more precise, this limit is reached in the weak topology of iJ 2 (ft). Indeed, it is a consequence of 
(|5.9[) that w^max € C„,([0, T,5 imax ]; ff 2 (ft)). Thus, — 1 + e < u(x) < 1 — e for all a; € ft and we can use u 
as a new "initial" datum and extend the solution (us limix , wg^max) beyond the time Ts tma , x . Moreover, 
the extension is still a classical solution since u is Holder continuous and uniformly separated from 
— 1 and 1. This contradicts the maximality of Ts, mlix and of (w^max, Wi5,max)- Hence, we necessarily 
have that Ts ima , x — T. 

To conclude the proof, we need to show that we can take the limit 5 \ and obtain a 
weak solution to Problem (P). With this purpose we rename simply as (usjWs) the maximal solution 
obtained in the previous part (which is now defined in the whole (0,T)), and observe that, thanks to 
estimates (|4.1I) . (I4.4j) . (|4.29[) . and (|4.22[) . there hold the following convergence relations: 



u s ^u weakly star in H 1 ^^-^') C\ L°°(0,T;V) C\ L°°((0,T) x ft), (5.15) 

eu s ->■ eu weakly in fl" x (0, T; ), (5.16) 

weakly in L 2 (0,T;V), (5.17) 

za^z weakly star in L°° (0,T;V) (1 L 2 (0,T; H 2 (n)), (5.18) 

w<5 = f(us) -> " weakly in i 2 (0, T; V), (5.19) 



for suitable limit functions u, w, z, v. The above properties, as well as the ones that will follow, are to 
be intended up to the extraction of (non- relabelled) subsequences of 5 \ 0. We then immediately see 
that relation (|l.ip passes to the limit. However, since it is only w € L 2 (0,T; V), the Laplace operator 
with the boundary condition d n w = have to be interpreted in the weak form through the operator 
A (cf. (12~T1) ). 

Next, applying the Aubin-Lions lemma, (|5.15[) gives 

u s ^u strongly in L p ((0,T) x ft), V>e[l,oo). (5.20) 

Thus, a standard monotonicity argument (see, e.g., [4j Prop. 1.1, p. 42]) permits to infer from (|5.19p 
that v = f{u) almost everywhere. Moreover, by the generalized Lebesgue's theorem, we get more 
precisely 

v s = f(u 5 ) -^v = f(u) strongly in i7((0,T) x ft), Vp e [1,2). (5.21) 
Now, as a consequence of estimate (I4.77j) . we obtain 

Wv s Vw weakly in L p ((0,T) x ft), for some p > 2. (5.22) 

Collecting (|5"^0|) and (1Q2"|) . we infer 

y|Vw 5 | 2 ^$ weakly in L p ((0,T) x ft), for some p > 1, (5.23) 

and for a suitable limit function $. 
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We can now write equation (|2.17[) for the solution (u$,ws) and see that, thanks to the above 
convergence relations, all terms pass to the limit. Actually, what we obtain for 6 \ is 

w = —Av + v + $ - Am + eu t , (5-24) 

at least in the distributional sense. To conclude the proof, we need to identify the function $. To this 
aim, we notice that, by (|5.15[) . (|5.17p . (I5.23[) , and a comparison of terms in (|2.29p . 

|| - Av 5 ||i>((o,T)xn) < c, for some p > 1. (5.25) 

Thus, we get that, also in the limit, —Av (the distributional Laplacean) lies in L p ((0, T) x ft) for some 
p > 1. More precisely, thanks to the no- flux condition and to elliptic regularity, we deduce 

v s -» v weakly in L p (0, T; W 2 ' p {n)), for some p > 1, (5.26) 

and d n v — on (0,T) x T in the sense of traces. Coupling (|5.21l) and (I5.26|) we obtain strong 
convergence of Vva by interpolation. Indeed, we can use (for example) the Gagliardo-Nirenberg 
inequality (see again [22]) i n the form 

||V(v« - v)\\ L2i/wm < c\\D 2 (v s - v^l^Wvs - v\\ l ^ 2/7[n) + || V - «|| il(Q) . (5.27) 

Then, integrating in time, and using the boundedness of W 2,p -norms resulting from (|5.26p . we infer 

||V(«<5 - t ')|| L i(0 : T;L 24 / 19 (O)) - C \\ D ( Vs ~ W )|Il 1 ((0 : T)xO)II W < 5 _ W ll/ 1 (0 : T;L 12 / 7 (O)) + H W<5 _ V II L!((0,T) xfi) 

<■ (wri 1 n 1 / 2 j-im 2 II 1 / 2 \\\ — n 1 / 2 I I _ II 

^ C\\\U "(5||£l((o,y)xf2) ' 11-^ W IIli((0,T)xO)/H 1; < 5 V W L 1 (0,T;L^f {Q.)) + W Vs W lli 1 ((0,T)xn) 

1/2 || 
< c\\v S - v \\ L i(Q^ T . L i2/7 {n)) + \\VS - w || L i((o,T)xn)- ( 5 ' 28 ) 

Owing to (|5.2ip . we then obtain that Vvs tends to Vu, say, strongly in L 1 ((0, T) x il), and consequently 
almost everywhere. Thus, recalling (|5.22p and using once more the generalized Lebesgue theorem, we 
have 

Vv 5 -> Vv strongly in L p ((0,T) x ft), for some p > 2. (5.29) 
whence, by (|5.20p . (|5.23p is improved up to 

y|Vu 5 | 2 ^|Vu| 2 strongly in L P ((0,T) x ft), for some p > 1. (5.30) 

Consequently, the function $ in (|5.24p is identified to its expected limit. Hence, we get (|2.29p (holding 
as a relation in L p (fl) for some p > 1, hence almost everywhere) and the boundary condition (|2.31|) . 
Finally, we notice that the limit 6 \ can be taken trivially in the initial condition (jl.3p . This 
concludes the proof. 

Remark 5.1. Theorem 12.41 states that uniqueness holds for weak solutions that are classical on all 
intervals (r, T), r > 0. This is, in fact, the same regularity class where we are able to prove existence. 
However, we cannot exclude that uniqueness might instead fail as one considers the (larger) class 
of all weak solutions, which in particular may contain some trajectory that does not achieve the 
"classical" regularity for strictly positive times. Actually, some of the calculations given in the proof 
(in particular, those related to the gradient terms) seem not be justified under the sole regularity 
conditions proper of weak solutions (which, for instance, may not be "separated" from the singular 
values ±1). 
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